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A discrete-event simulation approach is reviewed that does not require the knowledge of the solution of 
■ the wave equation of the whole system, yet reproduces the statistical distributions of wave theory by gen- 

^ I crating detection events one-by-one. The simulation approach is illustrated by applications to a two-beam 

interference experiment and two Bell test experiments, an Einstein-Podolsky-Rosen-Bohm experiment with 
single photons employing post-selection for pair identification and a single-neutron Bell test interferometry 
experiment with nearly 100% detection efficiency. 
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1 Introduction 

1 . 1 Quantum theory and the observation of single events 

Quantum theory has proven to be extraordinarily powerful for describing the statistical properties of a vast 
number of laboratory experiments. Starting from the axioms of quantum theory it is, at least conceptually, 
straightforward to calculate numbers that can be compared with experimental data as long as these numbers 
refer to statistical averages. However, if an experiment records individual clicks of a detector a fundamental 
problem appears. Although quantum theory provides a recipe to compute the frequencies for observing 
events it does not account for the observation of the individual events themselves, a manifestation of the 
quantum measurement problem [1,2]. A recent review of various approaches to the quantum measurement 
problem and an explanation of it within the statistical interpretation is given in Ref. [3] . From the viewpoint 
of quantum theory, the central issue is how it can be that experiments yield definite answers. As stated 
by Leggett [4], "In the final analysis, physics cannot forever refuse to give an account of how it is that 
we obtain definite results whenever we do a particular measurement". In fact, a common feature of all 
probabilistic models is that they do not entail a procedure that specifies how particular values of the random 
variables are being realized, but only contain a specification of the probabilities with which these values 
appear. As a probabilistic theory, quantum theory takes a special position in that it postulates that it is 
fundamentally impossible to go beyond the description in terms of probability distributions. At the present 
time, there is no scientific evidence that supports this assumption other than that it seems unsurmountable 
to account for the observation of the individual events within the context of quantum theory proper [1,2]. 
This suggests that the search for a cause-and-effect description of the observed phenomena should be done 
outside the realm of quantum theory. 
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1 .2 From single events to probability distributions and not vice versa 

In general, the problem can be posed as follows: Given a probability distribution of observing events, 
can we construct an algorithm which runs on a digital computer and produces events with frequencies that 
agree with the given distribution without the algorithm referring, in any way, to the probability distribution 
itself. An affirmative answer to this primary question can give rise to secondary questions such as "Can 
this be done efficiently?", "Is the algorithm unique?", and so on. 

There are several elements in the statement-of-the-problemthat deserve attention. The first is that we use 
"event" in the every-day meaning of the word. Thus, an event can be a click of a detector, the experimenter 
pushing a button, and so on. Second, a digital computer is nothing but a macroscopic physical device 
that evolves in time by changing its state in a well-defined, cause-and-effect, discrete-event manner, as 
specified by the algorithm. Although not practical, it is possible to build a mechanical apparatus that 
performs exactly the same function as the digital computer. Thus, one can view a simulation on a digital 
computer as a perfectly controlled experiment on a macroscopic mechanical system, simulating, in a cause- 
and-effect manner, the phenomenon that is being observed in the laboratory. Third, and most importantly, 
the algorithm should not rely in any way on the probability distribution of the events that it is supposed 
to generate. Otherwise it would be straightforward to use pseudo-random numbers and generate events 
according to this distribution. However, this is not a solution to the posed problem as it assumes that the 
probability distribution of the quantum mechanical problem is known, which is exactly the knowledge that 
we want to generate without making reference to quantum theory. Put differently, the algorithm should be 
capable of generating events according to an unknown probability distribution. 

In order to clarify our aim we draw an analogy [5] with the Metropolis Monte Carlo method for sim- 
ulating classical statistical mechanics [6], a primary example of a method that samples from an unknown 
probability distribution. According to the theory of equilibrium statistical mechanics, the probability that 
a system is in the state with label n is given by pn = e ''^'V S^=i e^^^" = e^^^"^ /Z, where N is the 
number of different states of the system, which usually is very large. En is the energy of the state, and 
(3 = l/ksT where ks is Boltzmann's constant and T is the temperature. Disregarding exceptional cases 
such as the two-dimensional Ising model, for a nontrivial many-body system the partition function Z is 
unknown. Hence, p„ is not known. We can now pose the question "Can we construct a simulation algo- 
rithm that generates states according to the unknown probability distribution (pi, . . . ,pjv)?" As already 
mentioned, an affirmative answer to this question was given a long time ago by Metropolis et al. [6]. The 
basic idea is to construct a dynamical system, a Markov chain or master equation that samples the space of 
N states such that in the long run, the frequency with which this system visits the state n approaches p„ 
with probability one. 

1.3 The event-based simulation approach 

The basic ideas of the simulation approach that we review in the present paper are that (i) we stick to 
what we know about the experiment, that is we consider the experimental configuration, its outcome and 
its data analysis procedure as input for constructing the simulation algorithm; (ii) we try to invent a set of 
simple rules that generates the same type of data as those recorded in an experiment while reproducing the 
averages predicted by quantum theory; (iii) we keep compatibility with macroscopic concepts. Our event- 
based simulation approach is unconventional in that it does not require knowledge of the wave amplitudes 
obtained by first solving a wave mechanical problem. Evidently, as outlined below, our event-based ap- 
proach requires a departure from the traditional way of describing physical phenomena, namely in terms 
of locally causal, modular, adaptive, classical (non-Hamiltonian) dynamical systems. 

Our event-based approach has successfully been used for discrete-event simulations of the single beam 
splitter and Mach-Zehnder interferometer experiment of Grangier et al. [7] (see Refs. [8-10]), Wheeler's 
delayed choice experiment of Jacques et aJ. [11] (see Refs. [10, 12, 13]), the quantum eraser experiment of 
Schwindt et al. [14] (see Ref. [10, 15]), two-beam single-photon interference experiments and the single- 
photon interference experiment with a Fresnel biprism of Jacques et al. [16] (see Ref. [10, 17]), quantum 
cryptography protocols (see Ref. [18]), the Hanbury Brown-Twiss experiment of Agafonov et al. [19] 
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(see Ref. [10,20]), universal quantum computation (see Ref. [21,22]), Einstein-Podolsky-Rosen-Bohm 
(EPRB)-type of experiments of Aspect et aJ. [23, 24] and Weihs et al [25] (see Refs. [5, 10, 26-30]), and 
the propagation of electromagnetic plane waves through homogeneous thin films and stratified media (see 
Ref. [10, 31]). An extensive review of the simulation method and its applications is given in Ref. [10]. In- 
teractive demonstration programs, including source codes, are available for download [32-34] . A computer 
program to simulate single-photon EPRB experiments can be found in Ref. [27]. So far the event-based 
simulation approach has not been used to simulate the phenomena of diffraction and evanescent waves. 

Given the above list of successful simulations, clearly an affirmative answer has been given to the 
primary question posed in section 1.2. A more detailed look to the algorithms that have been designed in- 
dicates that they are not unique. Regarding the efficiency of the algorithms we can say that (i) they are more 
efficient than their experimental counterparts because idealized models of the experimental equipment are 
considered and (ii) their efficiency is limited by the computational power of existing digital computers. Al- 
though the event-based simulation approach can be used to simulate a universal quantum computer [21,22], 
the so-called "quantum speed-up" cannot be obtained. This by itself is no surprise because the quantum 
speed-up is the result of a mathematical construct in which each unitary operation on the state of the quan- 
tum computer is counted as one operation and in which preparation and read-out of the quantum computer 
are excluded. Whether or not this mathematical construct is realized in Nature is an open question. 

1 .4 Event-based modeling and detection efficiency 

Applications in quantum information have increased the interest in single-particle detectors. High de- 
tection efficiencies are essential in for example quantum cryptography and some Bell test experiments. 
Single-particle detectors are often complex devices with diverse properties. In our event-based simulation 
approach we model the main characteristics of these devices by very simple rules. So far, we have de- 
signed two types of detectors, simple particle counters and adaptive threshold devices (see section 2.1). 
The adaptive threshold detector can be employed in the simulation of all single-photon experiments we 
have considered so far but is absolutely essential in the simulation of for example the two-beam single 
photon experiment. The efficiency, which is the ratio of detected to emitted particles, of our model detec- 
tors is measured in an experiment with one single point source emitting single particles that is placed far 
away from the detector. According to this definition, the simple particle counter has an efficiency of 100% 
and the adaptive threshold detector has an efficiency of nearly 100%. Hence, these detectors are highly 
idealized versions of real single photon detectors. No absorption effects, dead times, dark counts or other 
effects causing particle miscounts are simulated. 

In laboratory experiments, the single particle detection process is often quantified in terms of the ratio of 
detected to emitted particles, the overall detection efficiency. Evidently, the efficiency of the detector plays 
an important role in this overall detection efficiency but is not the only determining factor. Also the exper- 
imental configuration in which the detector is used plays an important role. Therefore, the experimenter 
usually choses the best overall performing single-particle detector for her or his particular experiment. Also 
in the event-based approach the experimental configuration plays an important role in the overall detection 
efficiency. Although the adaptive threshold detectors are ideal and have a detection efficiency of nearly 
100%, the overall detection efficiency can be much less than 100% depending on the experimental con- 
figuration. For example, using adaptive threshold detectors in a Mach-Zehnder interferometry experiment 
leads to an overall detection efficiency of nearly 100%, while using the same detectors in a single-photon 
two-beam experiment (see section 2.1) leads to an overall detection efficiency of about 15% [10, 17]. 

Also the data processing procedure which is applied after the data has been collected may have an 
influence on the final detection efficiency. For example, in the single-photon EPRB experiment of Weihs 
et al. a post-selection procedure with a time-coincidence window is employed to group photons, detected 
in two different stations, into pairs [25]. Clearly, this procedure omits detected photons and therefore 
reduces the final ratio of detected to emitted photons. This is also the case in the event-based simulation of 
this experiment (see section 2.2). Although simple particle counters with a 100% detection efficiency are 
used and thus all emitted photons are accounted for during the data collection process, the final detection 
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efficiency is less than 100% because some detection events are omitted in the post-selection data procedure 
using a time-coincidence window. 

In conclusion, even if ideal detectors with a detection efficiency of 100% would be commercially avail- 
able, then the overall detection efficiency in a single-particle experiment could still be much less than 
100% depending on (i) the experimental configuration in which the detectors are employed and (ii) the 
data analysis procedure that is used after all data has been collected. 

1 .5 Event-based modeling and the interpretations of quantum theory 

This paper is not about interpretations or extensions of quantum theory. The fact that there exist sim- 
ulation algorithms that reproduce the statistical results of quantum theory has no direct implications for 
the foundations of quantum theory. The average properties of the data may be in perfect agreement with 
quantum theory but the algorithms that generate such data are outside of the scope of what quantum the- 
ory can describe. Nevertheless, one could say that the event-based simulation approach is in line with 
the ensemble-statistical interpretation of quantum theory but also goes beyond this interpretation since the 
method is able to give a logical cause-and-effect description of how the ensemble is generated event by 
event. This makes metaphysical interpretations [35, 36] superfluous, at least for those experiments that can 
be simulated by the event-based approach. 

Of course, providing a description that goes beyond the statistical properties comes at a price. As it is 
well-known, "non-contextuality", literally meaning "being independent of the (experimental) measurement 
arrangement" is one of the properties that makes quantum (Maxwell's) theory so widely applicable [1]. To 
go beyond a statistical description unavoidably requires contextuality [1]. Therefore, the discrete-event 
approach provides a complete description of individual events but non-contextuality is lost. Note that as 
long as no measurement is performed the event-based description is non-contextual. For example, in the 
simulation of the single-photon EPRB experiment (see section 2.2) the two photons leaving the source 
have random but opposite polarization. The photons have a well-defined predetermined (non-contextual) 
polarization which corresponds to a certain vector S. When an observer makes a measurement on one 
photon using a polarizing beam splitter and two detectors, the observer gets a response +1 or -1 depend- 
ing on the angle of the polarizing beam splitter. Hence, depending on the angle of the polarizing beam 
splitter the observer measures +1 or -1 while the polarization vector S of the photon has a well-defined 
predetermined length and orientation. Thus, in this case, the contextuality (dependence on measurement 
arrangment) stems from the fact that the observer can only obtain partial information about the vector S 
when making a measurement. 

Finally, it should be noted that although the discrete-event algorithm can be given an interpretation as a 
realistic cause-and-effect description that is free of logical difficulties and reproduces the statistical results 
of quantum theory, at present the lack of relevant data makes it impossible to decide whether or not such 
algorithms are realized by Nature. Only new, dedicated experiments that probe more than just the statistical 
properties can teach us more about this intriguing question. Proposals for such experiments are discussed 
inRefs. [17,37]. 

2 Illustrative examples 

As a detailed discussion of an extensive set of discrete-event rules cannot be fitted in this short review, 
we have opted to present a detailed description of the algorithms for some fundamental experiments in 
quantum physics, namely interference of two coherent beams of particles and two Bell test experiments, 
an Einstein-Podolsky-Rosen-Bohm (EPRB) experiment with single photons [25] and a Bell's inequality 
interference experiment with single neutrons [38]. Our motivation to select the two-beam experiment with 
single photons as an example is that this experiment with minimal equipment shows interference in its 
purest form, that is without diffraction being involved. The experiment demonstrates that single particles 
coming from two coherent beams can gradually build up an interference pattern when the particles arrive 
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Fig. 1 Left: schematic diagram of a two-beam experiment with light sources So and Si of width a, separated by a 
center-to-center distance d. These sources emit coherent, monochromatic light according to the light flux distribution 
J{x, y) = S{x) [e(a/2 -\y - d/2\) + Q{a/2 -\y + d/2\)], where e(.) denotes the unit step function and with a 
uniform angular distribution, /3 denoting the angle. The light is recorded by detectors D positioned on a semi-circle 
with radius X and center (0, 0). The angular position of a detector is denoted by 6. Right: detector counts (circles) as 
a function of 9 as obtained from the event-based model simulation of the two-beam interference experiment depicted 
on the left. The solid line is a least-square fit of the simulation data to the prediction of wave theory, Eq. (5), with only 
one fitting parameter. Simulation parameters: on average, each of the 181 detectors receives lO'* particles, 7 = 0.99, 
a = c//, d = 5c//, X = 100c//, where c denotes the velocity and / the frequency of the particles. 



one by one at a detector screen. We choose the Bell test experiments as another example because the quan- 
tum theoretical description of this type of experiments involves entanglement. The single-photon EPRB 
experiment demonstrates that the two photons of a pair, post-selected by employing a time-coincidence 
window, can be in an entangled state. The neutron interferometry experiment shows that it is possible 
to create coiTelations between the spatial and spin degree of freedom of neutrons which, within quantum 
theory, cannot be described by a product state meaning that the spin- and phase-degree-of-freedom are en- 
tangled. In this experiment the neutrons are counted with a detector having a very high efficiency (« 99%), 
thereby not suffering from the so-called detection loophole. 

Furthermore, we have chosen to present the close-to-simplest algorithms that reproduce the results of 
quantum theory, hoping that this will help the reader grasp the basic ideas. For an extensive review of an 
event-based model for quantum optics experiments, see Ref. [10]. 

2.1 Two-beam interference 

In 1924, de Broglie introduced the idea that also matter, not just light, can exhibit wave-like properties [39]. 
This idea has been confirmed in various double-slit experiments with massive objects such as electrons [40- 
43], neutrons [44,45], atoms [46,47] and molecules such as Cqq and Cjo [48,49], all showing interference. 
In some of the double-slit experiments [ 1 6,4 1 ,42] the interference pattern is built up by recording individual 
clicks of the detectors. According to Feynman, the observation that the interference patterns are built up 
event-by-event is "impossible, absolutely impossible to explain in any classical way and has in it the heart 
of quantum mechanics" [50]. 

Reading "any classical way" as "any classical Hamiltonian mechanics way", Feynman's statement is 
difficult to dispute. However, taking a broader view by allowing for dynamical systems that are outside the 
realm of classical Hamiltonian dynamics, it becomes possible to model the gradual appearance of inter- 
ference patterns through a discrete-event simulation that does not make reference to concepts of quantum 
theory. As a discrete-event simulation is nothing but a sequence of instructions that changes the state of a 
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macroscopic classical apparatus (most conveniently a digital computer) in a prescribed manner, the demon- 
stration that we give in this section provides a rational, logically consistent, common-sense explanation of 
how the detection of individual objects that do not interact with each other can give rise to the interference 
patterns that are being observed. As there are several breeds of such models [10, 17, 51] that reproduce the 
statistical distribution predicted by quantum theory, only a new kind of experiment, specifically addressing 
this issue can provide a verdict about the applicability of these models to the problem at hand [17]. 

For the purpose of illustration, we consider the simple experiment sketched in Fig. l(left). This two- 
beam experiment can be viewed as a simplification of Young's double-slit experiment in which the slits 
are regarded as the virtual sources 5*0 and 5*1 [52]. In the two-beam experiment interference appears in 
its most pure form because in contrast to the two-slit experiment the phenomenon of diffraction is absent. 
The event-based model for this experiment has two types of events only: (i) the creation of one particle 
at one of the sources and (ii) the detection of that particle by one of the detectors forming the screen. We 
assume that all these detectors are identical and cannot communicate among each other and do not allow 
for direct communication between the particles, implying that this event-by-event model is locally causal 
by construction. Then, if it is indeed true that individual particles build up the interference pattern one by 
one, just looking at Fig. l(left) leads to the logically unescapable conclusion that the interference pattern 
can only be due to the internal operation of the detector [53]. Detectors that simply count the incoming 
photons are not sufficient to explain the appearance of an interference pattern and apart from the detectors 
there is nothing else that can cause the interference pattern to appear. Making use of the statistical property 
of quantum theory one could assume that if a detector is replaced by another one as soon as it detects one 
photon, one obtains a similar interference pattern if the detection events of all these different detectors are 
combined. However, since there is no experimental evidence that confirms this assumption and since our 
event-based approach is based on laboratory experimental setups and observations we do not consider this 
being a realistic option. Thus, logic dictates that a minimal event-based model for the two-beam experiment 
requires an algorithm for the detector that does a little more than just counting particles. 

2.1.1 Event-based model 

We now specify the model in sufficient detail such that the reader who is interested can reproduce our 
results (a Mathematica implementation of a slightly more sophisticated algorithm [17] can be downloaded 
from the Wolfram Demonstration Project web site [34]). 

• Source and particles: The particles leave the source one by one, at positions y drawn randomly from 
a uniform distribution over the interval [-d/2 - a/2, -d/2 + a/2] U [+d/2 - a/2, +d/2 + a/2]. The 
particle is regarded as a messenger, traveling in the direction given by the angle /3, being a uniform 
pseudo-random number between —tt/2 and tt/2. Each messenger carries with it a message e{t) = 
(cos27r/t,sin27r/t) that is represented by a harmonic oscillator which vibrates with frequency / 
(representing the "color" of the light). The internal oscillator is used as a clock to encode the time of 
flight t. When a messenger is created, its time of flight is set to zero. This pictorial model of a "photon" 
was used by Feynman to explain quantum electrodynamics [54] . The event-based approach goes one 
step further in that it specifies in detail, in terms of a mechanical procedure, how the "amplitudes" 
that appear in the quantum formalism get added together 

The time of flight of the particles depends on the source-detector distance. Here, we discuss as an 
example, the experimental setup with a semi-circular detection screen (see Fig. 1 (left)) but in principle 
any other geometry for the detection screen can be considered. The messenger leaving the source at 
(0, y) under an angle (3 will hit the detector screen of radius X at a position determined by the angle 
9 given by sin9 = {ycos^ j3 + sin (5^/ X'^ — j/^ cqs^ (5)/ X, where \y/X\ < 1. The time of flight is 
then given by t = ■\/ X'^ — 2yX sin 6 + y'^/c, where c is the velocity of the messenger The messages 
e{t) together with the explicit expression for the time of flight are the only input to the event-based 
algorithm. 
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• Detector: Microscopically, the detection of a particle involves very intricate dynamical processes [3]. 
In its simplest form, a light detector consists of a material that can be ionized by light. This signal is 
then amplified, usually electronically, or in the case of a photographic plate by chemical processes. 
In Maxwell's theory, the interaction between the incident electric field E and the material takes the 
form P • E, where P is the polarization vector of the material [52]. Assuming a linear response, 
P(a;) = x('^)E(a;) for a monochromatic wave with frequency uj, it is clear that in the time domain, 
this relation expresses the fact that the material retains some memory about the incident field, xi'^) 
representing the memory kernel that is characteristic for the material used. 

In line with the idea that an event-based approach should use the simplest rules possible, we reason as 
follows. In the event-based model, the fcth message — (cos 27r/tfe, sin 27r/tfc) is taken to represent 
the elementary unit of electric field E(<). Likewise, the electric polarization P(i) of the material is 
represented by the vector = {po,k,Pi,k)- Upon receipt of the kth message this vector is updated 
according to the rule 

Pfc =7Pfc-i + (1 -7)efc, (1) 

where < 7 < 1 and A: > 0. Obviously, if 7 > 0, a message processor that operates according to 
the update rule Eq. (1) has memory, as required by Maxwell's theory. It is not difficult to prove that 
as 7 — !• 1~, the internal vector converges to the average of the time-series {ei, 62, . . .} [10, 17]. 
The parameter 7 controls the precision with which the machine defined by Eq. (1) learns the average 
of the sequence of messages ei, 62, . . . and also controls the pace at which new messages affect the 
internal state of the machine [8]. Moreover, in the continuum limit (meaning many events per unit of 
time), the rule Eq. (1) translates into the constitutive equation of the Debye model of a dielectric [17], 
a model used in many applications of Maxwell's theory [55]. The learning process of the detector, 
governed by the rule Eq. (1) can namely be viewed as a process that proceeds in discrete time steps r, 
sothatpfc — p(A:r) — p(t) ande^ — e(fcr) = e{t). For small t, Pfc-i = p{t) — Tdp{t) / dt + 0{t^) 
and hence 

^--i^p(t) + ^eW. (2) 
ot r7 r7 

Letting the time step r, that is the time between the arrival of successive messages, approach zero 
(r — )• 0), letting 7 approach one (7 — > 1~) and demanding that the resulting continuum equation 
makes sense, leads to the following relation between r and 7, 

lim lim ^—^ = r, (3) 

T^OJ-^l- r7 

with < r < cxD. Hence, dp{t)/dt = —rp{t) +re{t) or in Fourier space P(a;) = x(w)E(a;), being 
the constitutive equation of the Debye model of a dielectric, with Eq. (3) giving an explicit expression 
for the relaxation time 1 /F in terms of the parameters of the event-based model. 

After updating the vector pfc, the processor uses the information stored in p^ to decide whether or not 
to generate a click. As a highly simplified model for the bistable character of the real photodetector 
or photographic plate, we let the machine generate a binary output signal 5*^ according to 

Sk = e{pl - n), (4) 

where 9(.) is the unit step function and < < 1 is a uniform pseudo-random number Note 
that the use of random numbers is convenient but not essential [10]. Since in experiment it cannot be 
known whether a photon has gone undetected, we discard the information about the Sk ^ detection 
events and define the total detector count as = ^j^ where k is the number of messages 

received. N is the number of clicks (one's) generated by the processor. 
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The efficiency of the detector model is determined by simulating an experiment that measures the de- 
tector efficiency, which for a single-photon detector is defined as the overall probability of registering 
a count if a photon arrives at the detector [56]. In such an experiment a point source emitting single 
particles is placed far away from a single detector. As all particles that reach the detector have the 
same time of flight (to a good approximation), all the particles that arrive at the detector will carry 
the same message which is encoding the time of flight. As a result pk (see Eq. (1)) rapidly converges 
to the vector corresponding to this identical message, so that the detector clicks every time a photon 
arrives. Thus, the detection efficiency, as defined for real detectors [56], for our detector model is 
very close to 100%. Hence, the model is a highly simplified and idealized version of a single-photon 
detector However, although the detection efficiency of the detector itself may be very close to 100%, 
the overall detection efficiency, which is the ratio of detected to emitted photons in the simulation of 
an experiment, can be much less than one. This ratio depends on the experimental setup. 

• Simulation procedure: Each of the detectors of the circular screen has a predefined spatial window 
within which it accepts messages. As a messenger hits a detector, this detector updates its internal 
state p, (the internal states of all other detectors do not change) using the message and then gener- 
ates the event Sk. In the case Sk ~ 1, the total count of the particular detector that was hit by the fcth 
messenger is incremented by one and the messenger itself is destroyed. Only after the messenger has 
been destroyed, the source is allowed to send a new messenger This rule ensures that the whole sim- 
ulation complies with Einstein's criterion of local causality. This process of creating and destroying 
messengers is repeated many times, building up the interference pattern event by event. 

2.1.2 Simulation results 

In Fig. 1 (right), we present simulation results for a representative case for which the analytical solution 
from wave theory is known. Namely, in the Fraunhofer regime {d <^ X), the analytical expression for the 
Ught intensity at the detector on a circular screen with radius X is given by [52] 



where ^ is a constant, q = 27r/ / c denotes the wavenumber with / and c being the frequency and velocity 
of the light, respectively, and 6 denotes the angular position of the detector D on the circular screen, see 
Fig. l(left). Note that Eq. (5) is only used for comparison with the simulation data and is by no means 
input to the model. From Fig. 1 (right) it is clear that the event-based model reproduces the results of wave 
theory and this without taking recourse of the solution of a wave equation. 

As the detection efficiency of the event-based detector model is very close to 100%, the interference 
patterns generated by the event-based model cannot be attributed to inefficient detectors. It is therefore 
of interest to take a look at the ratio of detected to emitted photons, the overall detection efficiency, and 
compare the detection counts, observed in the event-based model simulation of the two-beam interference 
experiment, with those observed in a real experiment with single photons [16]. In the simulation that 
yields the results of Fig. 1, each of the 181 detectors making up the detection area is hit on average by ten 
thousand photons and the total number of clicks generated by the detectors is 296444. Hence, the ratio of 
the total number of detected to emitted photons is of the order of 0.16, two orders of magnitude larger than 
the ratio 0.5 x lO^'^ observed in single-photon interference experiments [16]. 

2.1.3 What is the working principle? 

In our event-based approach the simple particle counters and the adaptive threshold detectors are ideal 
detectors with a detection efficiency of (nearly) 100%. In the case of the two-beam interference experiment 
considered in section 2. 1, using simple particle counters would not result in an interference pattern. These 
detectors simply produce a click for each incoming photon and do nothing with the information encoded 




(5) 
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in the messages e carried by the particles. These messages contain information about the time of flight 
of the particles, that is about the distance travelled by the particles from one of the two sources to one of 
the detectors constituting the circular detection screen. It is precisely the difference in the times of flight 
(or the phase differences) which is important in the generation of an interference pattern. Since, in the 
single-photon two-beam experiment the detectors are the only apparatuses available in the experiment that 
can process this information (there are no other apparatuses present except for the source) we necessarily 
need to employ an algorithm for the detector that exploits this information in order to produce the clicks 
that gradually build up the interference pattern. A collection of about two hundred independent adaptive 
threshold detectors defined by Eq. (1) and Eq. (4) and each with a detection efficiency of nearly 100% is 
capable of doing this. As pointed out earlier, the reason why, in this particular experiment, this is possible 
is that not every particle that impinges on the detector yields a click. 

Note that to simulate the interference pattern observed in single-photon Mach-Zehnder interferometry 
experiments it is possible, but by no means necessary, to use the adaptive threshold detectors which do 
not necessarily produce a click for each incoming photon [10]. Indeed, it suffices to use simple particle 
counters which produce a click for each incoming photon [8,9]. In the simulation of the Mach-Zehnder 
interferometry experiment the beam splitters process the information about the time of flight of the parti- 
cles [8-10]. Whether or not the detectors also process this information has no influence on the generation 
of the interference pattern [8-10]. Hence, in this case the number of detected photons is equal to the 
number of emitted photons in the simulation. 

2.2 Einstein-Podolsky-Rosen-Bohm (EPRB) experiment with single photons 

The EPRB experiment with photons, carried out by Weihs et al. [25, 57], is taken as a concrete example to 
illustrate how to construct an event-based model that reproduces the predictions of quantum theory for the 
single and two-particle averages for a quantum system of two spin- 1/2 particles in the singlet state and a 
product state [10, 30], without making reference to concepts of quantum theory. Recall that the quantum 
theoretical descriptions of the EPRB experiment with photons or with spin- 1/2 particles are identical. 

In short, the experiment goes as follows. A source emits pairs of photons. Each photon of a pair travels 
to an observation station in which it is manipulated and detected. The two stations are assumed to be 
identical. They are separated spatially and temporally, preventing the observation at station 1 (2) to have 
a causal effect on the data registered at station 2 (1) [25]. As the photon arrives at station i ~ 1,2, it 
passes through an electro-optic modulator (EOM) which rotates the polarization of the photon by an angle 
depending on the voltage applied to the modulator. These voltages are controlled by two independent 
binary random number generators. A polarizing beam splitter sends the photon to one of the two detectors. 
The station's clock assigns a time-tag to each generated signal. 

The firing of a detector is regarded as an event. At the nth event, the data recorded on a hard disk 
at station i ~ 1,2 consists of x„ ^ = ±1, specifying which of the two detectors fired, the time tag t„ ^ 
indicating the time at which a detector fired, and the two-dimensional unit vector a„ ^ that represents the 
rotation of the polarization by the EOM at the time of detection. Hence, the set of data collected at station 
i = 1,2 may be written as 



In the experiment, the data {Ti , T2} is analyzed after the data of a particular run has been collected [25]. 
Adopting the procedure employed by Weihs et al. [25, 57], we identify coincidences by comparing the 
time differences f„ 1 — tm,2 with a window W where n = 1, . . . ,Ni and m ~ 1, . . . , N2- By definition, 
for each pair of rotation angles a and b of the EOMs, the number of coincidences between detectors Dx,i 
(x = ±1) at station 1 and detectors Dy,2 (y = ±1) at station 2 is given by 



^ i {^71,2 tn,i^ ^n,i |^ — 1; • ■ • i • 



(6) 




(7) 



71—1 771^1 
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where 9(.) denotes the unit step function. In Eq. (7) the sum over all events has to be carried out such 
that each event (= one detected photon) contributes only once. Clearly, this constraint introduces some 
ambiguity in the counting procedure as there is a priori, no clear-cut criterion to decide which events at 
stations i ~ 1 and i ~ 2 should be paired. One obvious criterion might be to choose the pairs such that Cxy 
is maximum [57] but, such a criterion renders the data analysis procedure (not the data production) acausal. 
It is trivial though to analyse the data generated by the experiment of Weihs et aJ. such that conclusions do 
not suffer from this artifact [58]. The correlation of the two dichotomic variables x and y is defined as 

E^/ . \ + C — - C+- - C-+ 

where the denominator is the sum of all coincidences. In general, the two-particle averages E{a, b), and 
the total number of the coincidences not only depend on the directions a and b but also on the time window 
W used to identify the coincidences. For later use it is expedient to introduce the function [59] 



S = Sia, b, a', b') = E{a, b) - E{a, b') + E{a' , b) + E{a' , b'). (9) 

Local-realistic treatments of the EPRB experiment usually assume that the expression for the correla- 
tion, as measured in the experiment, is given by [60] 

N 
n=l 

which is obtained from Eq. (7) (in which each photon contributes only once) by assuming that N = Ni = 
N2, pairs are defined by n = m, and by taking the limit W ^ 00. However, the working hypothesis that 
the value of W should not matter because the time window only serves to identify pairs may not apply to 
real experiments. The analysis of the data of the experiment of Weihs et al. shows that the average time 
between pairs of photons is of the order of 30 /iS or more, much larger than the typical values (of the order 
of a few nanoseconds) of the time-window W used in the experiments [57]. In other words, in practice, the 
identification of photon pairs does not require the use of W's of the order of a few nanoseconds. The small 
value of W is required to maximize \S\ and is not really required for the data to violate Bell's inequality 
\S\ < 2. In other words, depending on the value of W, chosen by the experimenter when analyzing the 
data, the inequality 15*1 < 2 may or may not be violated. Hence, also the conclusion about the state of the 
system depends on the value of W, which turns W into a so-called context parameter Analysis of the data 
of the experiment by Weihs et al. shows that W can be as large as 150 ns for the Bell inequality to be 
violated [58] and in the time-stamping EPRB experiment of Agiiero et al. [61] \S\ < 2 is clearly violated 
for W < 9 fis. Hence, the use of a time-coincidence window does not create a "loophole". Nevertheless, 
very often it is mentioned that these single-photon Bell test experiments suffer from the fair sampling 
loophole, being the result of the usage of a time window W to filter out coincident photons or being the 
result of the usage of inefficient detectors [62]. The detection loophole was first closed in an experiment 
with two entangled trapped ions [63] and later in a single-neutron interferometry experiment [38] and in 
an experiment with two entangled qubits [64]. However, the latter three experiments are not Bell test 
experiments performed according to the CHSH protocol [59] because the two degrees of freedom are not 
manipulated and measured independently. 

The narrow time window W in the experiment by Weihs et al. mainly acts as a filter that selects pairs 
of which the individual photons differ in their time tags by the order of nanoseconds. The possibility that 
such a filtering mechanism can lead to correlations that are often thought to be a characteristic of quantum 
systems only was, to our knowledge, first pointed out by P. Pearle [65] and later by A. Fine [66], opening 
the route to a description in terms of locally causal, classical models. A concrete model of this kind was 
proposed by S. Pascazio who showed that his model approximately reproduces the correlation of the singlet 
state [67] with an accuracy that seems beyond what is experimentally achievable to date. Larson and Gill 
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showed that Bell-like inequalities need to be modified in the case that the coincidences are determined by 
a time-window filter [68], and models that exactly reproduce the results of quantum theory for the singlet 
and uncorrelated state were found [10,26,28,30]. Here, we closely follow Refs. [28,30]. 

2.2. 1 Event-based model 

A minimal, discrete-event simulation model of the EPRB experiment by Weihs et aJ. requires a spec- 
ification of the information carried by the particles, of the algorithm that simulates the source and the 
observation stations, and of the procedure to analyze the data. 

• Source and particles: Each time, the source emits two particles that carry a vector S„ ^ — (cos(^„ + 
(i — l)7r/2), sin(^„ + (* ^ 1)^/2)), representing the polarization of the photons. This polarization is 
completely characterized by the angle ^„ and the direction i ~ 1, 2 to which the particle moves. A 
uniform pseudo-random number generator is used to pick the angle < < 2tt. Clearly, the source 
emits two particles with a mutually orthogonal, hence correlated but otherwise random polarization. 
Note that for the simulation of this experiment it is not necessary that the particles carry information 
about the phase 27r/ti_„. In this case the time of flight „ is determined by the time-tag model (see 
below). 

• Electro-optic modulator: The EOM in station i = 1, 2 rotates the polarization of the incoming particle 
by an angle a^, that is its polarization angle becomes j = EOMi(^„ + (* ~ l)7''/2,ai) = £,n + 
{i — l)7r/2 — ai symbolically. Mimicking the experiment of Weihs et al. in which ai can take the 
values a, a' and a2 can take the values b, 6', we generate two binary uniform pseudo-random numbers 
Ai = 0, 1 and use them to choose the value of the angles ai, that is ai = a(l — Ai) + a'Ai and 
Q2 = b{l - A2) + h'A2. 

• Polarizing beam splitter: The simulation model for a polarizing beam splitter is defined by the rule 



where < r„ < 1 are uniform pseudo-random numbers. It is easy to see that for fixed ^ = ^j', this 
rule generates events such that the distribution of events complies with Malus law. 

• Time-tag model: As is well-known, as light passes through an EOM (which is essentially a tuneable 
wave plate), it experiences a retardation depending on its initial polarization and the rotation by the 
EOM. However, to our knowledge, time delays caused by retardation properties of waveplates, being 
components of various optical apparatuses, have not yet been explicitly measured for single photons. 
Therefore, in the case of single-particle experiments, we hypothesize that for each particle this delay 
is represented by the time tag [28, 30] 



that is, the time tag is distributed uniformly (0 < r'^ < 1 is a uniform pseudo-random number) over the 
interval [0, X{£,'n J]. For A(^^ J = Tq sin"' 2^^ ^ this time-tag model, in combination with the model 
of the polarizing beam splitter, rigorously reproduces the results of quantum theory of the EPRB 
experiments in the limit W ^ [28, 30]. We therefore adopt the expression A(^^j J = To sin"' 2^^ ^ 
leaving only To as an adjustable parameter. 

• Detector: The detectors are ideal particle counters, meaning that they produce a click for each in- 
coming particle. Hence, we assume that the detectors have 100% detection efficiency. Note that also 
adaptive threshold detectors can be used (see section 1.4) equally well [10]. 




(11) 



t. 



(12) 
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Fig. 2 Simulation results produced by the locally causal, discrete-event simulation model using the time-tag model 
Eq. (12). Left: 5* as a function of 9 with a — 9, a' — n /4 + 6, b = ir /8 and b' — 3n/8 for a time window W = 2 ns. 
The solid line connecting the circles is the result — 2\/2 cos 20 predicted by quantum theory. The dashed lines represent 
the maximum value for a quantum system of two spin- 1/2 particles in a separable (product) state (|S| — 2) and in a 
singlet state (|51 = 2\/2), respectively. Right: 15"! as a function of W foi a ^ 0, a' = 7r/4, b = 7r/8 and b' = Stt/S. 
For W ~ 200ns, the system changes from "quantum-like" (jSI > 2) to "classical like" (151 < 2). Simulation 
parameters: for each value of 6, the number of pairs generated is 3 x 10" (roughly the same as in experiment [25]) and 
To = 2000 ns, the adjustable parameter in the time-tag model Eq. (12). 



• Simulation procedure: The simulation algorithm generates the data sets T^, similar to the ones ob- 
tained in the experiment (see Eq. (6)). In the simulation, it is easy to generate the events such that 
Ni = N2. We analyze these data sets in exactly the same manner as the experimental data are 
analyzed, implying that we include the post-selection procedure to select photon pairs by a time- 
coincidence window W. The latter is crucial for our simulation method to give results that are very 
similar to those observed in a laboratory experiment. Although in the simulation the ratio of detected 
to emitted photons is equal to one, the final detection efficiency is reduced due to the time-coincidence 
post-selection procedure. 

This algorithm fully complies with Einstein's criterion of local causality on the ontological level. Once 
the particles leave the source, an action at observation station 1 (2) can, in no way, have a causal effect on 
the outcome of the measurement at observation station 2 (1). 

2.2.2 Simulation results 

In Fig. 2 we present some typical simulation results for the function S, showing that the event-based model 
reproduces the predictions of quantum theory for the singlet state. The single-particle averages Ei{a) and 
E2 (6) (data not shown) are zero up to the usual statistical fluctuations and do not show any statistically 
relevant dependence on b or a, respectively. Additional results, mimicking quantum systems in the singlet 
and product state, as well as a rigorous probabilistic treatment of the model can be found in [30]. For 
W = 50ns (data not shown), we find \S\ — 2.62 which compares very well with the values between 
2 and 2.57 extracted from different data sets produced by the experiment of Weihs et aJ. In all cases, the 
distribution of time-tag differences (data not shown) is sharply peaked and displays long tails, in qualitative 
agreement with experiment [57]. 

From Fig. 2(right), it follows that a violation of the Bell inequality jS*] < 2 depends on the choice of W, 
a parameter which is absent in the quantum theoretical description of the EPRB thought experiment. There 
are two limiting cases for which S become independent of W. If —5- 00, it is impossible to let a digital 
computer violate the inequality jS"! < 2 without abandoning the rules of Boolean logic or arithmetic [69]. 
For relatively small W (W < 150ns), the inequality |5| < 2 may be violated. When — 5- the 
discrete-event models which generate the same type of data as real EPRB experiments, reproduce exactly 
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the single- and two-spin averages of the singlet state and therefore also violate the inequality \S\ < 2. 
Obviously, as the discrete-event model does not rely on any concept of quantum theory, a violation of the 
inequality | | < 2 does not say anything about the "quantumness" of the system under observation [69-7 1 ] . 
Similarly, a violation of this inequality cannot say anything about locality and realism [69-72]. Clearly, 
the event-based model is contextual. The fact that the event-based model reproduces, for instance, the 
correlations of the singlet state without violating Einstein's local causality criterion suggests that the data 
{xn,i,Xn,2} generated by the event-based model cannot be represented by a single Kolmogorov probability 
space. This complies with the idea that contextual, non-Kolmogorov models can lead to violations of Bell's 
inequality without appealing to nonlocality or nonobjectivism [73, 74]. 

2.2.3 Why is Bell's inequality violated? 

In Ref. [30], we have presented a probabilistic description of our simulation model that (i) rigorously 
proves that for up to first order in W it exactly reproduces the single particle averages and the two-particle 
correlations of quantum theory for the system under consideration; (ii) illustrates how the presence of 
the time-window W introduces correlations that cannot be described by the original Bell-like "hidden- 
variable" models [60]. 

The time-coincidence post-selection procedure with the time-window W filters out the "coincident" 
photons based on the time-tags t„.i thereby reducing the final detection efficiency to less than 100%, 
although in the simulation a measurement always returns a H-1 or -1 for both photons in a pair (100% 
detection efficiency of the detectors). Hence, even in case of a perfect detection process the data set that is 
finally retained consists only of a subset of the entire ensemble of correlated photons that was emitted by 
the source. In other words, the use of a time-coincidence window destroys the "fair-sampling hypothesis". 

We briefly elaborate on point (ii) (see Ref. [30] for a more extensive discussion). Let us assume that there 
exists a probability P{xi,X2,ti,t2\ai,a2) to observe the data {xi,X2,ti,t2} conditional on {q;i,q;2}- 
The probability P{xi,X2,ti,t2\ai, 02) can always be expressed as an integral over the mutually exclusive 
events ^1, ^2, representing the polarization of the photons 

P{xi,X2,ti,t2\ai,a2) = / / P{xi,X2,ti,t2\ai,a2,^i,£.2)P{£.i,£.2\ai,a2)d£,id^,2- 

(13) 

According to Eq. (11) and Eq. (12), in the probabilistic version of our simulation model, for each event, (i) 
the values of xi, X2, ti, and t2 are mutually independent random variables, (ii) the values of xi and ti (X2 
and ^2) are independent of 02 and ^2 (cti and ^1), (iii) ^1 and ^2 are independent of ai or a2- With these 
assumptions Eq. (13) becomes 

P{xi,X2,ti,t2\ai,a2) = 



(ii) 



(i) 



(iii) 



J J P{xi,ti\ai,a2,^i,(2)Pix2,t2\ai,a2,^i,^2) 



xP(Ci,6|ai,a2)dCif^6 

r27T r2iT 



J J ^(^i'*il"i'^i)^(^2,i2|a2,6) 



47r2 

xP(Ci,6|ai,"2)dCit^6 

r2n r2iT 



-2 -P(a;i|ai,Ci)P(ti|ai,a)i"(a;2|«2,6)-P(i2|a2,6) 



47r2 

xP(Ci,6|ai:a2)dCi'^6 

^27r i'2-iT 



— P(a;i|ai,Ci)P(ti|«i,Ci)P(x2la2,6) 



/o Jo 

xP{t2\a2,^2)P{^u^2)d^id^2, 



(14) 
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which is the probabilistic description of our simulation model. The probabilistic model Eq. (14) is identical 
to the local-realist model used by Larsson and Gill in their derivation of a CHSH inequality with time- 
coincidence restriction [68]. 

According to our simulation model, the probability distributions that describe the polarizers are given 
by P{xi\ai, = [1 + Xi cos 2{ai — and those for the time-delays ti that are distributed randomly 

over the interval [0, A(^i + (« - 1)^/2 - cti)] are given by = d{ti)9{H^i + (i - l)7r/2 - 

oii) ~ti)/ + (j — l)7''/2 — OLi). In the experiment [25] and therefore also in our simulation model, the 
events are selected using a time window W that the experimenters try to make as small as possible [57]. 
Accounting for the time window, that is multiplying Eq. (14) by a step function and integrating over all ti 
and the expression for the probability for observing the event (xi, X2) reads 

/'27r /'27T" 

P(.Ti, X2|ai,a2) = / / P(a;i|ai,^i)P(a;2|a2,6)p(Ci,6|ai, "2)^6^6: (15) 

JQ Jo 

where the probabiUty density C2|ai, 02) is given by 

(16) 

The simple fact that ^21"!, 0^2) 7^ C2) brings the derivation of the original Bell (CHSH) inequal- 
ity to a halt. Indeed, in these derivations it is assumed that the probability distribution for ^1 and ^2 does 
not depend on the settings ai or 02 [2, 60]. Larsson and Gill have shown that due to the filtering by the 
time-coincidence window, Eq. (15) satisfies a modified CHSH inequality in which the upperbound of 2 is 
to be replaced by a number that depends on the ratio of the number of pairs satisfying the time coincidence 
criterion and the total number of pairs [68]. For finite W, this upperbound can be larger than 2. Hence, it 
is to be expected that Eq. (15) can violate the original CHSH inequality. 

By making explicit use of the time-tag model (see Eq. (12)) it can be shown that [30] (i) if we ignore 
the time-tag information (W > Tq), the two-particle probability takes the form of the hidden variable 
models considered by Bell [60], and we cannot reproduce the results of quantum theory [60], (ii) if we 
focus on the case W the single-particle averages are zero and the two-particle average E{ai, 02) = 
— cos(q;i — 02). 

In summary, although our simulation model and its probabilistic version Eq. (14) involve local pro- 
cesses only, the filtering of the detection events by means of the time-coincidence window W can produce 
correlations which violate Bell-type inequalities [66-68]. Moreover, for W ^ our classical, local and 
causal model can produce single-particle and two-particle averages which are the same as those of the 
singlet state in quantum theory. 

2.3 Violation of a Bell inequality in single-neutron interferometry 

The single-neutron interferometry experiment of Hasegawa et al. [38], demonstrating that the correlation 
between the spatial and spin degree of freedom of neutrons violates a Bell-CHSH inequality, is taken to 
illustrate how to construct an event-based model that reproduces this correlation by using detectors that 
count every neutron and without using any post-selection procedure. 

A schematic picture of the single-neutron interferometry experiment is shown in Fig. 3. Incident neu- 
trons are passing through a magnetic -prism polarizer (not shown) that produces two spatially separated 
beams of neutrons with their magnetic moments aligned parallel (spin up), respectively anti-parallel (spin 
down) with respect to the magnetic axis of the polarizer which is parallel to the guiding field B. The 
spin-up neutrons impinge on a silicon-perfect-crystal interferometer [45]. On leaving the first beam split- 
ter, neutrons may or may not experience refraction. A Mu-metal spin-turner changes the orientation of 
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Fig. 3 Diagram of the single-neutron interferometry experiment to test a Bell inequality violation (see also Fig. 1 in 
Ref. [38]). 



the magnetic moment from parallel to perpendicular to the guiding field B. The result of passing through 
this spin-turner is that the magnetic moment of the neutrons is rotated by 7r/2 (— 7r/2) about the y axis, 
depending on the path followed. Before the two paths join at the entrance plane of beam splitter BS3 a 
difference between the time of flights (corresponding to a phase in the wave mechanical description) along 
the two paths can be manipulated by a phase shifter. The neutrons which experience two refraction events 
when passing through the interferometer form the O-beam and are analyzed by sending them through a 
spin rotator and a Heusler spin analyzer. If necessary, to induce an extra spin rotation of tt, a spin flipper is 
placed between the interferometer and the spin rotator. The neutrons that are selected by the Heusler spin 
analyzer are counted with a neutron detector (not shown) that has a very high efficiency (« 99%). Note 
that neutrons which are not refracted by the central plate of the Si single crystal (beam splitters BSl and 
BS2) leave the interferometer without being detected. 

The single-neutron interferometry experiment yields the count rate N{a,x) for the spin-rotation angle a 
and the difference x of the phase shifts of the two different paths in the interferometer [38]. The correlation 
E{a,x) is defined by [38] 



^ ^ N{a, x) + N{a + ^,x + ^)- N{a + ^, x) ~ N{a, x + tt) 

N{a,x) + N{a + TT,x + '^)+ N{a + T:,x)+ N{a,x + '^y 

If quantum theory describes the experiment it is expected that for the neutrons detected in the O-beam, 

E{a, x) = cos(a + x) implying |5| = E{a, x) — E{a, x') + E{a' , x) + E{a' , x') > 2 for some values 
of a, a', X and x' so that the state of the neutron cannot be written as a product of the state of the spin and 
the phase. Experiments indeed show that |5| > 2 [38,75]. 



2.3.1 Event-based model 

A minimal, discrete event simulation model of the single-neutron interferometry experiment requires a 
specification of the information carried by the particles, of the algorithm that simulates the source and the 
interferometer components, and of the procedure to analyze the data. 

• Source and particles: The particles (neutrons) leave the source one by one and carry a message e = 
(e*'^i cos 61/2, e*^^ sin6'/2) where V» = ^nft + S, for i = 1,2. Here 6i - S2 and 6 specify the 
magnetic moment of the neutron and t specifies its time of flight where / is a frequency which is 
characteristic for a neutron that moves with a fixed velocity v. In the presence of a magnetic field 
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Fig. 4 Left: correlation E{a, x) between spin and path degree of freedom as obtained from an event-based simulation 
of the experiment depicted in Fig. 3. Solid surface: E{a, x) = cos(a + x) predicted by quantum theory; circles: 
simulation data. The lines connecting the markers are guides to the eye only. Model parameters: reflection percentage 
of BSO, . . . , BS3 is 20% and 7 = 0.99. For each pair (a, x), four times 10000 particles were used to determine the 
four counts N{a, x), N{a + tt, x + i"), N{a, X + ^) ™d N{a + tt, x + i"). Right: same as figure on the left but 
7 = 0.55. 



B = {Bx, By, Bz), the magnetic moment rotates about the direction of B according to the classical 
equation of motion. 

• Beam splitters BSO, . . ., BS3: The model for the beam splitter uses the learning process governed 
by the rule Eq. (1). Exploiting the similarity between the magnetic moment of the neutron and the 
polarization of a photon, we use the same model for the beam splitter as the one used in Ref. [10] 
for polarized photons. The only difference is that we assume that neutrons with spin up and spin 
down have the same reflection and transmission properties, while photons with horizontal and vertical 
polarization have different reflection and transmission properties [52]. 

• Mu metal spin turner: This component rotates the magnetic moment of a neutron that follows the 
H-beam (O-beam) by 7r/2 (— 7r/2) about the y axis. 

• Spin-rotator and spin-Hipper: The spin-rotator rotates the magnetic moment of a neutron by an angle 
a about the x axis. The spin flipper is a spin rotator with a = tt. 

• Spin analyzer: This component selects neutrons with spin up, after which they are counted by a 
detector. The model of this component projects the magnetic moment of the particle on the z axis and 
sends the particle to the detector if the projected value exceeds a pseudo-random number r. 

• Detector: The detectors simply count each incoming particle, meaning that we assume that the detec- 
tors have 100% detection efficiency. This is an idealization of the real neutron detectors which have a 
detector efficiency of more than 99% [76]. 



2.3.2 Simulation results 

In Fig. 4(left) we present simulation results for the correlation E{a, x), assuming that the experimental 
conditions are very close to ideal. For the ideal experiment quantum theory predicts that E{a,x) = 
cos(q: + x). As shown by the markers in Fig. 4, disregarding the small statistical fluctuations, there is 
close-to-perfect agreement between the event-based simulation data and quantum theory. The laboratory 
experiment suffers from unavoidable imperfections, leading to a reduction and distortion of the interference 
fringes [38]. In the event-based approach it is trivial to incorporate mechanisms for different sources of 
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imperfections by modifying or adding update rules. However, to reproduce the available data it is sufficient 
to use the parameter 7 to control the deviation from the quantum theoretical result. For instance, for 
7 = 0.55 the simulation (see Fig. 4) yields Smax = S{a = 0, x = 7r/4, a' = 7r/2, x' = 7r/4) = 2.05, in 
excellent agreement with the value 2. 052±0.010obtained in experiment [38]. For7 = 0.67 the simulation 
yields Smax — 2.30, in excellent agreement with the value 2.291 ± 0.008 obtained in a similar, more recent 
experiment [75]. 

2.3.3 Working principle 

From Ref. [10] we know that the event-based model for the beam splitter produces results corresponding 
to those of classical wave or quantum theory when applied in interferometry experiments. Important for 
this outcome is that the phase difference x between the two paths in the interferometer is constant for a 
relatively large number of incoming particles. If, for each incoming neutron, we pick the angle x randomly 
from the same set of predetermined values to produce Fig. 4, an event-based simulation with 7 = 0.99 
yields (within the usual statistical fluctuations) the correlation E{a, x) ~ [cos(a + x)]/2, which does not 
lead to a violation of the Bell-CHSH inequality (results not shown). Thus, if the neutron interferometry 
experiment could be repeated with random choices for the phase shifter x for each incident neutron, and the 
experimental results would show a significant violation of the Bell-CHSH inequality, then the event-based 
model that we have presented here would be ruled out. 

3 Discussion and outlook 

We have given a brief introduction to a new methodology for simulating, on the level of single events, what 
are usually considered to be "quantum" phenomena. In spirit, our approach is similar to cellular autonoma 
modeling advanced by Wolfram [77] or, for instance, lattice Boltzmann modeling of fluid dynamics [78]. 
The general idea is that simple rules, which are not necessarily derived from classical Hamiltonian dy- 
namics, define discrete-event processes which may lead to the (complicated) behavior that is observed in 
experiments. The basic strategy in designing these rules is to carefully examine the experimental procedure 
and to devise rules such that they produce the same kind of data as those recorded in experiment, while 
avoiding the trap of simulating thought experiments that are difficult to realize in the laboratory. The event- 
based model is entirely classical in the sense that it uses concepts of the macroscopic world and makes no 
reference to quantum theory but is nonclassical in the sense that some of the rules are not those of classical 
Newtonian dynamics. 

Depending on the experimental configuration, the simulation of the interference and correlation phe- 
nomena, which are observed when individual photons and neutrons are detected one by one, can make use 
of one or two of the four following detection processes: 

• (1) The detectors are simple particle counters producing a click for each incoming particle and no 
post-selection data procedure using a time-coincidence window is used. The final detection efficiency 
is 100%. 

• (2) The detectors are simple particle counters and a post-selection data procedure using a time- 
coincidence window is used. Although the ratio of detected to emitted particles is one, the final 
detection efficiency is less than 100%. 

• (3) The detectors are adaptive threshold detectors not producing a click for each incoming particle 
and no post-selection data procedure using a time-coincidence window is used. The final detection 
efficiency can be less than 100%, depending on the experimental configuration. 

• (4) The detectors are adaptive threshold detectors and a post-selection data procedure using a time- 
coincidence window is used. The final detection efficiency is less than 100%. 
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Table 1 Overview of tlie event-based simulation of plioton (P) and neutron (N) experiments classified according 
to the detection process: (1) Detectors are simple particle counters and no post-selection data procedure using a 
time-coincidence window; (2) detectors are simple particle counters and post-selection data procedure using a time- 
coincidence window; (3) detectors are adaptive threshold detectors and no post-selection data procedure using a time- 
coincidence window; (4) detectors are adaptive threshold detectors and post-selection data procedure using a time- 
coincidence window. Numbers in brackets are references to the papers in which simulation results are reported. The 
symbol "V" ("X") indicates that it is possible (impossible) to obtain the results of quantum theory but that we did not 
publish the simulation results. The symbol "NA" denotes that the particular detection process is not applicable to the 
specified experiment. 



Experiment 


Particle 


(1) 


(2) 


(3) 


(4) 


Single beam splitter experiment 


P 


[8,9] 


NA 


[10] 


NA 


Mach-Zehnder interferometry experiment 


P 


[8,9] 


NA 


[10] 


NA 


Wheeler's delayed choice experiment 


P 


[12,13] 


NA 


[10] 


NA 


Quantum cryptography protocols 


P 


[18] 


NA 


V 


NA 


Quantum eraser experiment 


P 


[15] 


NA 


[10] 


NA 


Single photon tunneling 


P 


V 


NA 


[10] 


NA 


Two-beam interference experiment 


P 


X 


NA 


[10,17] 


NA 


Reflection and refraction from an interface 


P 


[31] 


NA 


[10] 


NA 


Multiple beam fringes with a plane-parallel plate 


P 


[31] 


NA 


[10] 


NA 


Quantum computation 


P 


[21,22] 


NA 


V 


NA 


Einstein-Podolsky-Rosen-Bohm experiment 


P 


X 


[5,26-30] 


X 


[10] 


Classical correlations in HBT experiment 


P 


X 


X 


[20] 


X 


Quantum correlations in HBT experiment 


P 


X 


V 


X 


[10] 


Interferometry 


N 


Section 2.3 


NA 


V 


NA 


Bell-type experiment; path-spin correlations 


N 


Section 2.3 


NA 


V 


NA 


Interferometry; stochastic and deterministic absorption 


N 


V 


NA 


V 


NA 


Interferometry; path-spin-energy entanglement 


N 


V 


NA 


V 


NA 


Interferometry; time-dependent blocking of one path 


N 


V 


NA 


V 


NA 



In Table 1 we give an overview of which type of experiments our event-based approach can simulate with a 
particular detection process. From Table 1 it can be clearly seen that discarding detection events in a post- 
selection data procedure using a time-coincidence window or using detectors which do not produce a click 
for each incoming particle is not a characteristic of the event-based simulation approach for simulating 
interference and quantum correlation phenomena. 

Three fundamental experiments, two with photons and one with neutrons, were used to illustrate our 
approach. On purpose, we have kept the event-based model as simple as possible, perhaps creating the 
impression that each experiment will require its own set of rules. This is not the case. One universal 
event-based model for the interaction of photons with matter suffices to explain, without altering the rules, 
the interference and correlation phenomena that are observed when individual photons are detected one 
by one [10]. This universal model produces the frequency distributions for observing many photons that 
are in full agreement with the predictions of Maxwell's theory and quantum theory [10]. The same model 
was used for simulating a single neutron interferometry experiment. Similarly, we have constructed an 
event-based model that simulates a universal quantum computer, reproducing the probability distributions 
of the quantum mechanical system without actually knowing them [22]. 

An important question is whether an event-based model leads to new predictions that may be tested 
experimentally. In the stationary state (after processing many events) the event-based model reproduces 
the statistical distributions of quantum theory. Therefore, new predictions can only appear when the event- 
based model is operating in the transient regime, before the event-based model reaches its stationary state. 
Experiments with a Mach-Zehnder interferometer consisting of two independent 50% reflective beam split- 
ters and a phase shifter (not an integrated Mach-Zehnder interferometer) and two-beam interference that 
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may be able to address this issue have been discussed in Ref. [37] and Ref. [17], respectively. We hope that 
our simulation work will stimulate the design of new experiments to test the applicability of our approach 
to event-based processes and to exclude some of our non-unique models. 

Finally, it may be of interest to mention that the discrete-event approach reviewed in this paper may 
open a route to rigorously include the effects of interference in ray-tracing software. For this purpose, it is 
necessary to extend the event-based model to include diffraction, evanescent waves and the behavior of the 
Lorentz model for the response of material to the electromagnetic field [55]. We leave these extensions for 
future research. 
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